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Abstract 

We present two classes of exact solutions to a geometric model which describes the 
kinetics of fragmentation of d-dimensional hypercuboid-shaped objects. The first class 
of exact solutions is described by a fragmentation rate a{xi, ...,Xd) = 1 and daughter 

distribution function b{xi, ..,Xd\x[, ■■■,x'^) = ■ ■ . . The second class of 

exact solutions is described by a fragmentation rate a(xi, x^) = xi"^ ...Xrf"'*/2'^ and 
a daughter distribution function ..,Xd\x'i, ■■■■,x'^) = 2'^6{xi — x'i/2). ..5{xd — x'j2). 
Each class of exact solutions is analyzed in detail for the presence of scaling solutions 
and the occurrence of shattering transitions; the results of these analyses are also 
presented. 



PACS: 0520, 0250. 



1 Introduction 

Fragmentation occurs in numerous important physical (droplet break-up [1] and fiber length 
reduction [2]), chemical (depolymerization through, shear action [3,4], chemical attack [5] 
and exposure to nuclear, ultra-violet and ultra-sonic radiation [6,7]) and geological (rock 
crushing and grinding (communition) [2]) processes. Theoretical predictions of the evolution 
with time of the size distributions of the fragmenting objects during such processes is of 
great interest and importance. There are essentially two approaches in use for determining 
the evolution in time of the object size distributions as a function of the initial conditions 
and the fragmentation rates. The first approach relies upon statistical and combinatorial 
arguments [8-10]. The second approach has been through the analysis of the kinetic equation 
modelling the fragmentation [11-13]. 

In the kinetic equation approach the fragmentation process can be described by the 
evohition in time of the size distribution c{x, t), where x is the size of the fragments and t is 
the time, through a kinetic equation. This theoretical approach is of a mean field nature since 
fluctuations are ignored. Fragments are assumed to be distributed homogeneously at all times 
throughout the system, i.e., there is perfect mixing and the shape of the fragments is ignored. 
Consequently, the size of the fragments is the only dynamical variable that characterizes a 
fragment in the kinetic equation approach. A number of aiithors have expended much effort 
in finding exact solutions to the kinetic equation, in order to study specific practical problems 
and to provide a greater understanding of the behaviour of physical, chemical and geological 
systems in which fragmentation occurs [5,14-21]. Although the basic kinetic equations are 
linear, and in principle soluble, the number of exact solutions is few, mainly because of the 
non-local structure of the kinetic equations. 

Of considerable importance are the scaling solutions. These are essentially the solutions 
in the long-time {t — > oo), small-size {x — > 0) limit where the size distribution c{x,t) evolves 
to a simpler universal form. This form is universal in the sense that it does not depend on 
the initial conditions. Most experimental systems evolve to the point where this behaviour 
is reached. Scaling theories based on a linear kinetic equation have been derived for a large 
class of models which undergo fragmentation [18, 22-24]. 

The time evolution of the fragmentation process depends qualitatively on the behaviour 
of the probability of the break-up for the fragments. For break-up rates increasing sufficiently 
quickly with decreasing size or mass, a cascading break-up occurs in which a finite part of 
the total size or mass is transfered to fragments of zero or infinitesimal size or mass. This 
so-called "shattering" [17, 25] or "disintegration" [26] phenomenon is accompanied by a 
violation of the usual dynamical scaling as well as violation of size or mass conservation. 

The general form of the 1-dimensional multiple fragmentation equation is given by 



where a{x) gives the rate of fragmentation of particles of size x, the daughter distribution 
function b{x\x') is the average number of particles of size x produced when a particle of size 
x' breaks up, and c{x,t) is the size distribution of particles of size x at time t. To ensure 
that the total size or mass of the fragmenting particles is conserved we insist that 
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= r dx'x'b(x'\x) (2) 

^0 



holds. The average number of particles produced when a particle of size x fragments is given 

by 

< N(x) >= r dx'b{x'\x). (3) 

^0 

On physical grounds we must have < N{x) >> 2, which together with (2) and (3), places 
constraints on the possible choices for b{x\x') available to us. 

In the special case of binary fragmentation where 2 particles are produced per fragmenta- 
tion event (1-3) can be rewritten in terms of the single symmetric function F{x, x') = F{x', x) 
as follows. Firstly, write, 

a(x) = r dx'F{x - x', x'). (4) 

Then to ensure that there are precisely 2 particles produced per fragmentation event, choose 

u/ I '\ '2F(^x,x — x) , , 

b{xx)^ — , 5 

a{x') 

in which case (1) becomes 

= -c(x, t) r dx'F(x - x', x') +2 r dx'F(x, x' - x)c(x\ t) (6) 

Ot Jo Jx 

where F{x,x') describes the rate at which particles of size {x + x') fragment into particles 
of size X and x'. 

As mentioned above, the kinetics of such 1-dimensional fragmentation processes is now 
well- understood with numcroTis explicit exact solutions, scaling solutions and quantitative 
descriptions of shattering transitions known to us. 

In reahstic fragmentation processes particles have both size and shape, and it is clear 
that this geometry of these fragmenting particles will influence the fragmentation process. 
For example a particle may be selected for fragmentation at a rate which is dependent on its 
area or volume, but the manner in which the fragmentation of the particle is implemented 
will, in general, depend on its precise dimensions. If 2 particles have the same area but 
one is needle-shaped and the other is square-shaped, they may be equally likely to fragment 
as a consequence of their possesing the same area, but the needle-shaped particle is much 
more likely to fragment across its longer side, whereas the square-shaped particle is equally 
likely to fragment across either side. Until recently, all these properties were represented by 
a single parameter, namely the size or mass of the fragmenting particle. 

Recently, various authors [27-29] have introduced and investigated simple kinetic models 
describing the fragmentation of 2-dimensional and more generally d-dimensional particles. 
These authors present several simple classes of explicit exact solutions for 2-dimensional 
models. Krapivsky and Ben-Naim [28] also discuss the presence of scaling and multi-scaling 
in their models of fragmentation for rf-dimensional particles. In [29], the shattering transition 
in a 2-variable fragmentation model is investigated. 

However, as far as we know there do not exist any explicit exact solutions for fragmenting 
particles in d-dimensions, for general d. We present in this paper, 2 classes of exact solutions 
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to a geometric model which describes the kinetics of fragmentation of dimensional objects, 
for general d. Each class of exact solutions is analyzed in detail for the presence of scaling, 
multi-scaling, and the occurrence of shattering transitions. 



2 Fragmentation in (/-dimensions 

The general form of the d-dimensional multiple fragmentation equation is given by 
dc(xi,...,Xd,t) 

— = -a{xi,...,Xd)c{xi,...,Xd,t) 

at 

noo rco 

+ / dx[... dx'da{x[,...,x'd)b{xi,...,Xd\x[,...,x'd)c{x[,...,x'd,t) (7) 

where c(xi, Xd^ t) is the size distribution of particles with size Xi...Xd at time t, a{xi^ ■■■,Xd) 
is the rate at which a particle characterized by {xi,...,Xd) fragments, and the daughter 
distribution function b{xi, ...,Xd\x[, ...,Xd) is the rate at which a particle characterized by 
{xi, ...,Xd) is produced from a particle characterized by {x[, 

Consequently, the average number of particles produced per fragmentation event is 

< N{xi,...,Xd) >= dx[... dx'db{x[,...,x'd\xi,...,Xd) (8) 
Jo JO 

and mass conservation requires, 

r^i ^ rxd ^ ^ 

Xi...Xd= dx^... dx'dx[...x'db{x[, ...,Xd\xi, ...,Xd). (9) 
Jo Jo 

Equations (7-9) form a complete set of equations, which define the fragmentation process 
given a{xi, ...,Xd), b{xi, ...,Xd\x[, ■■■,x'^) and suitable initial conditions. Of course, the func- 
tions a{xi, ...,Xd) and6(a;i, ...,Xd\x'i, ■■.,x'd) must be chosen to make the equations physically 
meaningful and furthermore b{xi, ...,Xd\x[, ■■■,x'd) must be chosen to ensure that (8) and (9) 
hold, as well as the obvious physical constraint N{xi, ...,Xd) > 2. 

As a special case, consider fragmentation processes such that a given rf-dimensional par- 
ticle fragments into 2^^ pieces per fragmentation event. In this case, as in section 1 above, we 
may rewrite (7-9) in terms of a single function F{xi,x[] ..;Xd,x'd) as follows. Firstly, choose 



rxi rXd 

a{xi, Xd) = / dx\... / dx'^F{xi - x\,x'^] ..; Xd - x'^, x'^) (10) 
Jo JO 



fXl l-Xd 

dx^... / 

lo Jo 

where F{xi,x[; ..;Xd,Xd) is the rate of fragmentation of a particle characterized by {xi + 
x[)...{xd + x'd) into 2^^ smaller particles characterized by: xi...Xd,..., x[...Xd. Now choose, 

7/ I / /\ Xi — Xi] Xd, x^ — Xd) 
b{xi,...,Xd\xi,...,Xa) ^ (11) 

a[Xi, Xd) 

where, of course, 

F{xi, Xi — xi] ..; Xk, Xf.; Xd, Xd — Xd) — F[xi, Xi — xi; ..; Xj^, Xk', Xd, Xd — Xd) (12) 
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i.e., F{xi, x[ — xi, ..; Xk, x'l^; Xd, x'^ — Xd) is symmetric in all pairs of arguments {xk, x'f.) for 
k = 1, ...,d. 

It is easily shown that < N{xi, ...,Xd) >= 2*^, as required, and that mass conservation 
requires, 

Xi-.-Xd — f^i H-r' H-r' Fl'-r _ ^' ^' ■ ■ ^ , _ ^' ' ^ ^ 

Jo '-l-^l-'-Jo UX^r l^Xi Xi,Xi,..,Xd X^,X^) 

The fragmentation equation (7) now becomes, 
dc{xi,...,Xd,t) /-^i 



dt 



-c(xi,...,Xd,t) dx[... dx'aF(xi- x[,x[;..;xd- x'a,x'a) 
Jo Jo 

roc roo 

+ 2'^ dx[... dx'aF{x[- xi,xi;...;x'a- Xd,Xd)c{x[,...,x'a,t) (U) 

J x\ Jx^ 



where F{xi,x\] ..]Xd-iX'^ is defined by (10-13). 
Model 1 

This model is described by the fragmentation rate and daughter distribution function given, 
respectively, by 

a{xi,...,Xd) = 1 (15) 

Ux. xAx' x') - K + 2)-K + 2)a:r-^? .^g. 
0[xi,...,Xd\Xi,...,Xa) - . 

Insisting that mass is conserved per fragmentation event places the following restrictions on 
the homogeneity indices cci, a^, 

ai>-2 ;i^l,...,d. (17) 

It is easily shown that the number of fragments produced per fragmentation event is 

( (ai+2)...(arf+2) . if all rv- > -1 

<Nix^,...,Xd)>^{ ' "aiia,> | (ig) 

[ oo ; II some or all ai < —1 

On physical grounds, for all «j > —1, we must insist that when < N(xi, ...,Xd) > is finite it 
must be so that < A^(a;i, ...,Xd) >> 2. This constrains the by 

^^i±%4^>2. (19) 
The d-dimensional multiple fragmentation equation (7) now becomes 

dc{xi,...,Xd,t) 

— ^ -c{Xi,...,Xd,t) 



dt 

+ ia^ + 2)...{ad + 2)xr...xT / -^c{x[, ...,x'„t).{20) 

Jxi Xi JXct X^ 
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We need to solve (20) subject to appropriate initial conditions, 

c{xi, Xd, 0) = /(xi, Xd) ^ 0. (21) 
Define the Laplace transform of c{xi, ...,Xd,t) with respect to t, 4>{xi, ...,X(i,s), by 

(t){xi,...,Xd,s) = I dte~^^c{x-i,...,Xd,t) (22) 
Jo 

in which case we may recover c{xi, ...,Xd,t) from the inverse Laplace transform 

]_ /■7+ioo 

c{xi,...,Xd,t) ^ — dse^^(t){xi,...,Xd,s) (23) 

ZTTZ J-y-ioo 

where 3fJ(s) > 7 to ensure convergence. 

Taking the Laplace transform of (20) with respect to t yields 



0(xi, ...,Xd,s) 
+ 



^ fixi,...,Xd) 

is + l) 

{ai + 2)...{ad + 2)x'^\..x'^'' dx[ 



s + 1) 



tol+l ••• 
xi Xi 



Xd"^ 



<l>(x[,...,x'„s).{2A) 



Following the approach introduced in [30] one finds 



(f){Xi, ...,Xd,s) 



+ 



^ f{xi,...,Xd) 

(s + l) 

{ai + 2)...{ad + 2)x'^\..x^'' dx[ 



is + l) 



00 -1 

X E 



=0 {r\y 



ai + 2)...{ad + 2) ^^^xU 

(s + 1) ^Xi^'" ^Xa 



(25) 



Performing a simple contour integration yields 



c{xi,...,Xd,t) =e H f {xi,..., Xd) + {ai + 2)... {ad + 2)x"K..x'^'' 



X 



X 



dx[ 

'ai+l ■■■ 



dx'. 



ai + 2)...(a, + 2)lnp)...lnp 

^Xi^ ^Xd 



,=-^(r + l)!(r!)'' 
For mono-disperse initial conditions, 

/(xi, Xd) = 5{xi - li)...5{xd - Id 

(26) becomes 



(26) 



(27) 
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c(xi,...,Xd,t) =e ^6{xi-li)...6{xd-ld) + {ai + '2)--{o:d + 2) 



x'^K..x 



/ai+l ]ad+^ 



X 



E 



f-r+l 



^(r + l)!(r!r 



(ai + 2)...K + 2)ln(^)...ln(^) 



(28) 



For ctj = and d = 1 (28) reduces to the result of Ziff and McGrady [16]. For ctj = 
and d — 2 (28) reduces to the solution presented by Rodgers and Hassan [27]. When the 
ctj are completely general and d = 1, (28) is equivalent to the exact solution for the model 
investigated by McGrady and Ziff [17], with P — —1. 



Model 2 



In this model we investigate a fragmentation rate and daughter distribution function 
given, respectively, by 

a{xi,...,Xd) = ^x^\..Xd'' (29) 

b{x„...,Xd\x[,...,x',) = 2%x^-^)...5{xd-^). (30) 

It is easily shown that mass conservation per single fragmentation event holds, and that the 
average number of particles per fragmentation event, < N{xi, ...,Xd) >, is 2*^. 

This particular choice for the fragmentation rate a{xi, Xd) and the daughter distribu- 
tion function b{xi, ...,Xd\x[, ■■■,Xd) can be implemented by F{xi,x[; ...;xd,x^) with 

Fixi,x[; ...■,Xd,x'd) = {xi + x[)°'\..{xd + XdT''S{xi - x[)...d{xd - x'd). (31) 

In this form the kinetics of the model become a little more transparent. The fragmentation 
rate F{xi,x[; ...-jXdjX'^) describe a fragmentation process in which an object splits into 2'^ 
equally (?) sized pieces. 

This choice for a{xi, Xd) and b{xi, Xd\x[, x'd), or equivalently, F{xi,x[; Xd, x'd), 
reduces the d-dimensional multiple fragmentation equation (7) to the following form, 

^"^"^^g-;"''^'^^ = -^^^c(xi, Xd, t) + 2'^-xT...xTc{2x„ 2xd, t) (32) 

where a — ai + ... + ad. 

Solving (32) subject to the initial conditions 

c{xi, Xd, 0) = f{xi, Xd) 7^ (33) 
via the approach outlined in [30], one finds 
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c{xi,...,Xd,t) = e ^i'-^^'*/^'|^/(a;i,...,a;d) 



me/,- 

with a = ai + ... + cifi 7^ and — {0, 1, 2, A; — 1, /c + 1, r}. 
Using the fact that 

i^oS n (2-«-2^") ) H ^^^^ 

we can obtain c{xi, Xd, t) for a — ai + ... + aa — without any extra effort. Exphcitly, 



cix„ ...,Xd,t) = e--i --/^/^^ Y: S^(xr-^^^)7(2'-xi, ...,rxd) (36) 

r=0 

where, of course, a — ai + ... + a^,. 
For mono-disperse initial conditions, 

f{xi, Xd) = 6{xi - li)...d{xd - Id) (37) 

we find 



c{xi, ...,Xd,t) = < 



e-'i'-'^'*/2'' 5(a;i - h)...5{xd - Id) + E 2^'+'-)^'''+''' 



r=l 



x5(xi — -^)...5{xd — ^) - — Y[ (2'"«_2fcc«) — 



A:=0 



r=0 



;q;^0 . (38) 
;q; = 



For d = 1, these results reduce to those presented in [30]. For d = 1 and a = ai = 0, one 
recovers the exact solution of Bak and Bak [19]. 



3 Scaling and multi-scaling 



To investigate the presence of scaling or multi-scaling, we introduce the d-tuple Mellin trans- 
form of the distribution function c{xi, ...,Xd,t) defined by 

roo roc 

M{si,...,Sd,t) = dxi... dxdxl'-\..x'/-^c{xi,...,Xd,t). (39) 
Jo Jo 

The functions M{si, Sd, t) for fixed Si, Sd are known as the moments of the distribution 
function c{xi, ...,Xd,t). 

Combining (7) and (39) gives 
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poo roo 

/ dxi... / dxda{xi,...,Xd)c{xi,...,Xd,t) 
Jo Jo 



dt 

x{x{^-\..xr' - I dx[...J^ dx'A''~'---^d' K<^---,^'d\^i,---,^d))- (40) 

We now investigate the the two classes of exact solutions to the d-dimensional multiple 
fragmentation equation, presented in section 2, for the presence of scaling and multi-scaling. 



Model 1 



In this case (40) becomes, 

dM(si,...,Sd,t) ^ / {ai + 2)...{ad + 2) 
dt \ {si + ai)...{sd + ad) 

provided Sj + ccj > for i = 1, d. Then it follows that 

M{si,...,Sd,t)^M{si,...,Sd,0)e T^T+^^Trj^l+^ (42) 

provided Sj + ccj > for i = 1, d. 

To obtain the number of objects, A^(i), in our fragmenting system we must take Si — 1 
for i = 1, d in (42). We find that 

/ (ai+2)...(arf+2) \ 

iV(t) = M(l,...,l,t) = M(l,...,l,0)e ^ (ai+i)...(ad+i)J* (43) 

which is only valid for a, > — 1 and i = l,...,d. When some or all of the ctj < —1, then 
(41) and (42) are not valid, since it can easily be shown that the number of particles in the 
system is infinite, even though the mass or volume, V{t), of the system is conserved. If we 
set Sj = 2 for i = 1, d in (42) we find that the volume, V{t), of our system is constant, 
i.e., 

\/(t) = M(2,...,2,t) = M(2,...,2,0) (44) 

provided ojj > —1 for i = I, d. 

An interesting feature of (41), for Sj+aj > and i — 1, d, is that it implies the existence 
of an infinite number of conservation laws, apart from the usual volume conservation one 
usually encountered. The moments M{si, Sd, t) with si, Sd satisfying 

+ f ^ + = 1 (45) 

{si + ai)...{sd + ad) 

are all time- independent. Of course, Si ^ 1 \/i — 1, ...,d in this case, otherwise (45) would 
contradict (19) even if ctj > —1 — l,...,d. Besides, the number of objects, N{t), in a 
fragmenting system cannot possibly be conserved. Thus, in addition to the conservation of 
the total volume, V{t), there are an infinite number of hidden conserved integrals for those 
Si, with i = l,...,d, that lie on the hypersurface defined by (45). According to Krapivsky 



\M{si,...,Sd,t) 



(41) 
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and Ben-Nairn [28], it is precisely these integrals which are responsible for the absence of 
scaling solutions for the d-dimensional fragmentation processes. Indeed, the scaling solution 

c(a;i, Xd, t) ~ r<P{fxi, ...,fxd) (46) 
implies an infinite number of scaling relations 

w — zs (47) 

where s — si + ... + s^, which together with (45) cannot all be satisfied by the scahng 
exponents w and z. This rules out the possibility of scaling solutions in this model, however 
multi-scaling solutions may be possible. 



Model 2 



For this case (40) becomes, 

dM{si,...,Sd,t) ^ _ 1^ 

dt [ 2('-<^) 2d 



\M{si + ai,...,Sd + ad,t) (48) 



where s = Si -|- ... -|- sa- Again, as in model 1, we observe that (48) implies the existence 
of an infinite number of conservation laws. The moments M{si, Sd,t) with the Sj, for 
i — 1, ...,d, satisfying 

s = 2d (49) 

where s = Si + ... + Sd, are all time-independent. Thus, in addition to the conservation 
of the total volume, V{t) = M(2, ...,2,t), there are an infinite number of hidden conserved 
integrals for those Sj, with i — 1, d, that lie on the hypersurface defined by (49). Due to 
the existence of an infinite number of hidden conserved integrals, we do not expect there to 
be any scaling solutions to this model. However, this docs not appear to be the case, and 
we demonstrate this explicitly by finding a scaling solution to this model. 
In this case assume a scaling solution of the form 

c(xi, Xd, t) ~ VityUixi/Vit), Xd/Vit)) (50) 

as i — > oo, where V{t) is the volume of a typical time-dependent cluster volume. This scaling 
form is fully equivalent to (46), however its form is more convenient for our investigations [18]. 
The exponent —2d ensures conservation of volume of the complete system of fragmenting 
particles. 

A short calculation will show that 

M(si, Sd, t) ~ Vity-^^mis,, Sd) (51) 
where the scaling moments m{si, Sd) are defined by 

poo roc 

m{s„...,Sd)^ dCi... / dUi''-'-4/'-'HCi,-,Cd). (52) 
Jo Jo 

If the moments are to be conserved, then they must be constant. This occurs when 
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s = 2d (53) 

which is in complete agreement with (49), as it should be. 
Substituting (50) into (32) yields 



dV{t) 

— t 

( - ^0(6, Cd) + 2^+x2ei, 2e.))e^^..c 



V{tf^' dt 



ad 
d 



(54) 



where the separation constant u is positive since V{t) must be a decreasing function of time 
in a fragmenting system and = Xi/V{t) for i — 1, ...,d. Then 



V{t) ~ < 



a > 0, i — > oo 

a — 0, i — > oo . (55) 
a < 0, t <tc 



[ (ic-i)H 

These expressions are only valid provided scaling holds. For a < a singularity is encoun- 
tered within a finite time t^, and scaling becomes invalid. In this instance, we anticipate 
shattering which will be discussed later. When a = 0, as in the 1-dimensional case [24, 26, 
30], a scaling form does not exist which is consistent with (f){^i, ^d) — > constant as — >■ 
and 0(^1, ^d) ^ as — > oo for A; = 1, d. Consequently, we need only look at the case 
when q; > 0, for which scaling is valid. 

When a > 0, we assume that 0(^i, ^d) vanishes at = and = oo for A; = 1, d. 
It can be shown from (54) that 

as ^fe — > oo for A; = 1, d. Therefore, we will assume that 

'^(^1' ^<^) = 772— —72^/(^1' -' ^<^) (57) 
1^1 + ••• + ^d) 



where we insist that /(^i, ^rf) = 1 at = oo for /c = 1, d. 

Sustituting (57) into (54) and performing a lengthy calculation yields a solution for 
/(6, -,Cit)- Hence, for a > 0, 

Ui + ••• + 4dj v n=i n (2"^" - 1) / 

m=l 

Thus, a rather unusual situation occurs in this d-dimensional model. A scaling solution 
exists, in spite of the fact that we have an infinite number of hidden conserved integrals, 
of the form (50), as t — > oo with V{t) given by (55) and (f){^i, ■■■,^d) given by (58). These 
results arc consistent with the exact solution of model 2 presented in section 2 above. When 
d = 1 these results reduce to those in [30] and are of a similar form to those of Cheng and 
Redner [24]. 



10 



4 Shattering transitions 



Formally, the volume, V{t), of the system is defined by 

roc roc 

V{t) ^ dxi... dxdXi...Xdc(xi,...,Xd,t), (59) 
Jo Jo 

so that , with the aid of (7) and (9) it can easily be shown that 

IM.o (ao) 

indicating that the volume, V{t), is conserved. However, when the fragmentation rate in- 
creases sufficiently fast as the volume of the fragments decreases to zero, a cascading of the 
fragmentation occurs such that volume is lost to fragments of zero or inifinitesimal volume. 
This cascading process, which has been named "shattering" [17, 25] or "disintegration" [26], 
is somewhat similar to gelation in coagulating systems, where mass is lost to an infinite gel 
molecule [31, 32]. Gelation and shattering are both signalled by the condition dV{t)/dt < 0. 
When shattering is suspected a more subtle analysis to that used to derive (60) is required. 
To analyze the shattering transition define a cut-off volume Ve{t), with < £ <^ 1, by 

roc roo 

Ve{t) = j dxi... j dxdXi...Xdc{xi,...,Xd,t), (61) 

with 

V{t) = \imV,{t). (62) 
It is easily shown that the cut-off volume loss is given by 



dV.it) 



dt 



dxi... J dxda{xi, Xd)c{xi, Xd, t) 

X / dx[... [ dx'^x[...x'ab{x[,...,x'a\xi,...,Xd) (63) 
Jo Jo 

where, of course, 

= ,ta (64) 

dt 6^0+ dt ^ ' 



Model 1 



In this model shattering does not occur for any values of the homogeneity indices cc^, for 
i = l,...,d. This is easily demonstrated from the definition of volume, V(t), of our system 
via (59) and the exact solution for this particular model (26). One finds that the volume, 
V{t), of our system is both finite and time- independent. 
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Model 2 



As can be easily demonstrated by substituting the exact solution for this particular model 
(36) into the definition of the volume, V{t), of the system (59), shattering does not occur 
when the sum of the homogeneity indices, a, is zero. Again, one finds that the volume, V{t), 
of the system is both finite and time-independent. It remains to investigate the case when 
the sum of the homogeneity indices, a, is less than zero. 

To analze the shattering regime, o; < 0, we need to know the behaviour of c{xi, ...,Xd,t) 
as i — > oo and — > such that txk remains fixed for k — 1, ...,d. A short calculation will 
show that the asymptotic form of c{xi, ...,Xd,t) in the shattering regime is 

c{xi,...,Xd,t)r.T{t)x^K..x^'' (65) 
where T{t) ^ 0. The exponents Ai, A,^ are restricted by 

A=|a|-2rf (66) 

where A = Ai + ... + A^. In this analysis it does not matter what each individual value of 
the Afe, for A; = 1, is, all that matters is what the sum. A, is. Therefore, we need not 
concern ourselves with explicitly determining the A^, for A; = 1, d. 
For q; < 0, (63) becomes 

^^ = -^y^ dxi...^ dx,xr+^..x^''+'c(xi,...,a:,,t) (67) 
for this particular model. Substituting (65) into (67) yields 

dl4(t) _ 1^,,, 
di 2'^ 

V ctj + Aj = —2, with \ <i <d and 



r(i)(ln2)'' (68) 



(69) 



dt 2'^^^^ U «^^^ + Afc+i + 2 ^'■■■V «^ + Arf + 2 

if ccj + Aj = —2 for i = 1, k and Oj + Aj 7^ — 2 for i = /c + 1, d, where < A; < d — 1. 
Hence, we see that 



di £-^0+ di 
which proves that shattering does indeed take place for o; < 0. 



d^W-U^lM)^0 (70) 



5 Volume distributions 



In this section we briefly consider the volume distribution function C(V, t) deflned by 

/•oo /"OO 

C{V,t)^ dxi... dxdS{xi...Xd-V)c{xi,...,Xd,t) (71) 
Jo Jo 

which is usually very useful for providing a partial description of a fragmenting system. 
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Model 1 



For = 0, i = 1, ...,d, and mono-disperse initial conditions the relevant exact solution 
to this model is 



c{xi,...,Xd,t) = e ^^6{xi - li)...6{xd 
1 



E 



ld) + 



(72) 



Substituting this into (71) yields 



C{V,t) 



ody. oo 



Kr+1)-1] 



(73) 



This reduces to the result of Ziff and McGrady [16] when d = 1, and to the result of Rodgers 
and Hassan [27] for d = 2. Expanding (73) for small t gives 



C{V,t) 



(74) 



(<i-l)! 

As d increases by one, the power of the logarithmic divergence in the second term in (74) 
also increases by a factor of one. Analogous to the 1-dimensional case, the d-dimensional 
case forms the borderline case for the shattering transition. 
Define the normalized nth moments of the volume, V{t), by 



/o°°d\/\/"C(\/,t) 

jrdvc{v,t) • 

Then, it follows that 

< >Vn^ g-2=*(l-l/(n+l)'')t/n 

which indicates that C{V,t) does not exhibit a scahng form. 



(75) 



(76) 



Model 2 

In this model, for mono-disperse initial conditions, the exact solution is given by (38). 
Substituting (38) into (71) yields 



C{V, t) - 



r=l k=0 11 ^ ' 



r=0 



;q; 7^ 
;a — 



(77) 
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These results may be compared with the results when d = 1 for c{xi, ...,Xd,t) given in (38) 
or in [30]. The similarity is quite remarkable, indicating that the scaling and shattering 
behaviour of C{V,t) given in (77) will match closely that observed for c{xi, ...,X(i,t) given 
by (38), when d—1. 

6 Conclusions 

In reality, fragmenting particles will have both size and shape, i.e., a geometry. Intrigued 
by the possibility that the geometry of the fragmenting particles may influence the frag- 
mentation process, we have investigated two distinct d-dimensional fragmentation mod- 
els, for d > 1. Two classes of exact solutions to these geometric models, which describe 
the kinetics of fragmentation of ci-dimcnsional particles, arc presented. The first class 
is described by a fragmentation rate a{xi, ...,Xd) = 1 and a daughter distribution func- 
tion b{xi,...,Xd\x[,...,x'a) = (ai + 2)...(ad + 2)a;i"i...a;d"Va;i'("i+^)...a;/"'^+^). For d > 1, 
this particular class of exact solutions exhibits multi-scaling and does not permit the oc- 
currence of a shattering transition. The second class of exact solutions is described by 
a fragmentation rate a{xi, ...,Xd) = xi°'^ ...Xd°''^ and a daughter distribution function 
b{xi, ..,Xd\x[, ...,x'^) = 2'^6{xi — x[/2)...6{xd — x'^/2). This particular class of exact solu- 
tions describes a type of fragmentation process in which particles always break-up into 2'^ 
equally sized pieces at various rates which depend upon the geometry of the fragmenting par- 
ticles and the homogeneity indices ai, ...,ad- this type of fragmentation has been observed 
and studied when polymers degrade under tension (stretching) [33], or in the presence of 
a destructive force-field such as ultra-sound [6]. Defining, a, to be the sum of all the ho- 
mogeneity indices ai, ad, it is shown that this particular class of exact solutions exhibits 
scaling for a > and this scaling form is explicitly determined. For a = 0, we show that 
multi-scaling is likely and that the shattering transition is not permitted. When o; < 0, we 
show that a shattering transition occurs. 

An interesting scenario occurs in our investigation into the scaling behaviour of the second 
class of exact solutions. When a > 0, we have shown that a scaling solution to our model 
exists and we explicitly find this scaling solution. This is very surprising, since it has been 
suggested [28] that scaling solutions are not supposed to exist when we have an infinite 
number of hidden conserved integrals. We suggest that the existence of an infinite number 
of hidden conserved integrals is not always indicative of the absence of scaling solutions. 
There must exist other important criteria which conclusively indicate the absence of scaling 
in a particular d-dimensional fragmention process. We propose to investigate the nature of 
these criteria in subsequent work. 

An investigation into the occurrence of a shattering transition in the second class of exact 
solutions presented in section 2 is quite intriguing, when a < it is shown that shattering 
occurs. It is interesting that whether shattering occurs or not is determined exclusively by 
the fact that a < 0, and not on the sign of the individual cti, ad, which add up to give a. 
Of course, the ferocity of the shattering transition will depend on how negative a is, with 
the possibility of competing effects between positive and negative ai, i — 1, ...,d, which will 
act to moderate the ferocity of this shattering transition. 

The volume distribution function, C{V, t), is very useful for providing a partial description 
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of a fragmenting system. For small t the 2-dimensional model differs from the 1-dimensional 
model by the presence of a logarithmic divergence term. The 3-dimensional model differs 
from the 1-dimensional model by the presence of a logarithm squared divergence term, and 
so on. An analysis of the normalized nth momentsof the volume V indicate that C{V, t) does 
not exhibit scahng. For the second class of exact solutions presented in section 2, it is shown 
that the similarity between C{V,t) and c{xi, ...,Xd,t), with d = 1, is quite remarkable. This 
is a strong indicator that the scaling and shattering behaviour of C{V, t) will be very similar 
to that observed for c{xi, ...,Xd,t), with d= 1. 

We have found that the introduction of more than one parameter to characterize the 
geometry of the fragmenting particle can have a significant effect on the kinetics of the frag- 
mentation process. As a propect for future research, one could investigate problems with 
sources and sinks in rf- dimensions, given that one now has a useful method available for de- 
termining exact solutions to geometric models which describes the kinetics of fragmentation 
of d-dimensional particles. 
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